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Fiber Coupler 
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We simplify the coupling theory between two contiguous, parallel, 
multimode step-index fibers, describe the coupling concept, and derive 
an upper estimate for the overall coupling efficiency between the two 
fibers. The maximum coupling obtainable, according to this estimate, 
is less than 72 percent (-1.5 dB). The coupling efficiency derived for 
short coupling lengths shows good agreement with experimental re- 
sults. 

I. INTRODUCTION 

A multimode fiber tap -coupler is a useful component for certain op- 
tical communication systems such as the optical data bus. Unfortunately, 
it is not an easy matter to evaluate the simultaneous coupling process 
between the hundreds or thousands of modes. An analysis of the problem 
has been given by Snyder and others. 1 " 3 These authors analyzed only 
the coupling between certain mode pairs. Snyder 4 recently reported the 
total power transition between multimode fibers. However, his conclu- 
sion is based on HEi m modes for analyzing the crosstalk. 

We have derived a similar simplified expression for the total coupling 
between identical, contiguous, parallel, step-index, multimode fibers, 
which can expand to all modes under the restriction that two fiber cores 
are touching each other. We predict that the maximum coupling effi- 
ciency is less than 1.5 dB when all modes carry an equal amount of power. 
The distance between two fibers affects the coupling efficiency very 
seriously when fibers have large numerical apertures. 

Our simple formula agrees very well with experimental results, in spite 
of a large number of approximations made. 

II. COUPLING COEFFICIENT 

2. 1 Simplified coupling coefficient 

Figure 1 shows the geometry of the fiber coupler. The cores are parallel 
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Fig. 1 — Geometry used to derive the coupling coefficients. 

to each other and surrounded by a medium that has the same refractive 
index as the cladding. 

A general description of strongly coupled multimode fibers is given 
in Appendixes A and B. The coupling coefficient in eq. (24) has a com- 
plicated form. It involves the modified Bessel functions and many dif- 
ferent parameters that depend upon the eigenvalue equations. We will 
simplify the coefficient under the following assumptions: (i) the two 
fibers are identical, (ii) only coupling between modes having the same 
propagation constant is considered (see discussion in Section 2.2), and 
(Hi) the distance, d, between the two fiber axes is nearly equal to 2a, 
where a is the core radius of the fiber. 

We can rewrite the coupling coefficient by using assumptions (i) and 
(ii). The result is 



'AiBi 



= \C 



BiAi 



V2R u 2 \K 2 i(wd/a) cos (2la) ± K (wd/a)\ 
7 a v 3 K[-i(w)Ki+i(w) 



where 

7 = l k when / = 
7 = 1 when / > 
a = (defined in Fig. 1) 
a = radius of core 
no = refractive index of core 
A = relative refractive index difference 
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d = distance between the two fibers 

/ = azimuthal order number 

u 2 = w 2 + u 2 

v 2 = w 2 + u 2 

u 2 = a 2 k 2 nl(2A). 

The coupling coefficient (1) is still complicated because it includes the 
modified Bessel functions. Further simplification is necessary to present 
a simple physical picture for the coupling process. 

We use assumption {Hi) to simplify (1) using the asymptotic expan- 
sions of the modified Bessel function. As discussed in Appendix C, the 
modified Bessel function term in (1) can be expressed by a very simple 
expression for a fairly large range of azimuthal order numbers. We use 
the average value of | C, | as shown in Appendix C. The result is as fol- 
lows: 

where 

d/a s2 

7 = 1 for all modes except LPo m ; y = l k for LPo m . 

We let 7 equal 1 for all modes. This assumption does not seriously 
affect the results of this analysis. 
The coupling coefficient is thus, 

\/2A /u 2 V2w 
where 



^I.Vp/KfvJS.,-/.-*,), (8) 

a \u d Vird/a I 



d/a s 2. 

The coupling coefficient expressed by (3) does not explicitly depend upon 
the azimuthal order number /. However, it is still dependent upon the 
solution of the eigenvalue equation. We simplify eq. (3) further by in- 
troducing simple expressions for u, w, and v. 

The parameter u is expressed by the total number of modes in the 
fiber when the total number of modes is large. The results are 5,6 

v 2 = 2N, (4) 

where 

v » 1 

N = total mode number of fiber A and B. 
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We order all LP modes according to their z -propagation constants 
from the largest to the smallest. We label them by the sequential num- 
bers, i = l-N. For example, the two orthogonal LP i modes are desig- 
nated as the first and the second mode. The LPn modes assume the 
orders 3, 4, 5, and 6. For i » 1, the cutoff value of u for the ith mode is 
approximately 5 

Cutoff =(20 1/2 (i»D- (5) 

We replace u by u cu toff and we obtain 

u = (2i) 1/2 (i » 1). (6) 

Equation (3) can be simplified by (4) and (6). The result is, 

ra-isi -£&(£) (i-i)" 

x exp[-(2N - 2i)** (dla - 2)] (?) 
Vd/a 

or 

2 3/4( A) i/4 , i v . i v 1/4 exp[-(2iV - 2t) 1/2 (d/a - 2)] 



m- 



~VTT.k-n a 3/2 \N/ \ Nt Vdfc 

where 

a = radius of core 
k = 2tt/X 

no = refractive index of core 
A = relative refractive index difference 
d = distance between the two fibers. 
When dla = 2, the coupling coefficient becomes, 

91/4 Al/4 /i\ / / \ 1/4 

^ = W=M^(n)( 1 -») ■ <8) 

This simple expression for the coupling coefficient does not require 
the eigenvalue solutions. Figure 2 shows that the coupling coefficient 
reaches a maximum for the mode number i = (4/5)N. The reason is the 
following: the coupling coefficient as expressed by (24) is based on the 
field interaction of the evanescent field tail of a mode of fiber A and the 
core field of the same mode order in fiber B. Generally speaking, the 
higher-order modes have a stronger field in the cladding relative to their 
fields in the core. Therefore, the field interaction between the field tail 
of a mode of fiber A and the core field of the same order in fiber B in (16) 
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Fig. 2 — Coupling coefficient vs mode order. 

can be expected to increase with mode order until it reaches a maximum. 
Very high mode orders have very small fields in the core and therefore 
the interaction between very-high-order modes decreases again. 

Figure 2 shows the effect of the gap between two fiber cores. Coupling 
occurs only between the higher-order modes as the gap increases. We 
will discuss the coupling efficiency using (7) and (8) in the next sec- 
tions. 

2.2 Coupling efficiency 

The power coupled from the ith mode of fiber A to the ith mode of 
fiber B can be obtained from (25) when the ith mode of fiber A carries 
unit power and the ith mode of fiber B carries zero power at z = 0. The 
coupled power is 



where 



{P A ^B)i = *in 2 [\Ci\z], 



z = the coupling length. 



(9) 



MULTIMODE FIBER COUPLER 733 



The total coupled power is obtained by the summation of (9). (The 
coupling coefficients for some of the modes are nearly equal to zero; 
especially equal to zero are the coupling coefficients for one of the modes 
with / = 0. So the coupling efficiency defined here is an upper bound.) 
If all modes of fiber A carry equal power, then the total power from fiber 
A to fiber B is, 

^zh Z[sinH\Ci\z)]. (10) 

Pin Ni=l 

If we treat the mode number i as a continuous variable, the coupling 
efficiency rj becomes 

Pa-B 



v = 



Pi, 



<f r C" [am* (\Ci\z)di\. (ID 

N Jo 

The coupling between the near synchronous modes is not negligible, 
especially when d /2a = 1. However, the upper bound of the coupling 
efficiency is expressed by (1 1) when | Cy 1 1 C n | « |C; 1 2 for i - a £ j ^ 
i + a (a/N) « 1, a defined by (25). When the two fiber cores are touching 
each other, the coupling efficiency r] is 

<N _ 

sin 2 \Ci\z di 



-£"&&*-*>* ,ra 



/ ( 2 1/4 A 1/4) Z 

= I sin' 
Jo 

where 



i 



Figures 3 and 4 show several examples. 

III. DISCUSSION 

As a first example, consider two fibers having the following parame- 
ters: NA = 0.2, a = 25 Mm (radius), n = 1.457, y = 0.633 urn. This fiber 
carries about 1,250 modes (N = 1,250). 

The coupling coefficient is 

— - = /_! \/na \ (0(1 ~ 1/4 exp \-V2N(VT=t)(d/a - 2)] 

According to Fig. 2, the maximum coupling efficiency (1.5 dB) is achieved 
at d/2a = 1 and 0.011 z/a = 3.67. Thus, the coupling length z is about 
334a or z = 8.36 mm. 
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Fig. 3 — Coupling efficiency vs coupling length. 
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Fig. 4 — Coupling efficiency vs coupling length. 
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Fig. 5— Coupling efficiency vs coupling length; points indicate experimental results. 

When d /2a = 1.01 so that the gap between the two fibers is about 0.5 
ixm, the coupling efficiency drops to 3.2 dB with the same coupling length 
(about 8.36 mm). When d/2a = 1.05 so that the gap between the two fi- 
bers is about 2.5 nm, the coupling efficiency becomes about 13.7 dB with 
the same coupling length. Thus, the distance between two fibers affects 
the coupling efficiency very strongly. 

We now look at another example with respect to the following fiber 
parameters: NA = 0.4, a = 50 ^m, no = 1.457, X = 0.633 /im. This fiber 
carries about 20,000 modes. The coupling coefficient is expressed by 

IC-I = 0.01- (0(1 -t) l/4 



x exp [-2 >/2N(VT=t)(d/2a - 1)] 

Vd/2a 

The maximum coupling efficiency (about 1.5 dB) is achieved when d/2a 
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SIMPLIFIED COUPLING COEFFICIENT 
LIMITATION OF COUPLING COEFFICIENT 
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Fig. 6 — Coupling coefficient vs mode order; points indicate numerical results. 

= 1 and the coupling length is 16.7 mm. However, if the ratio d/2a 
changes by only 1 percent so that the gap between the fibers becomes 
about 1 i*m, the coupling efficiency becomes about 11.7 dB. The gap 
between two fibers affects the coupling efficiency more seriously when 
the fiber has a large numerical aperture and large radius. 

Practically, it is not easy to produce a fiber coupler that has a uniform 
interaction gap over a long coupling length. Figure 5 shows experimental 
results that only cover very short coupling lengths. Plastic-clad fiber is 
attached to the acrylic base and its cladding is peeled off over the re- 
quired coupling length. The exposed cores are pushed together and form 
a parallel coupling region. Silicone of the same type as the fiber cladding 
is injected to form a common cladding around the coupling region. Figure 
4 shows that this simple theoretical approach yields good agreement. 
Figures 6 and 7 show results obtained from a numerical calculation of 
the coupling coefficient (23) and for fibers with the following parameters: 
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Fig. 7— Coupling efficiency vs coupling length; points indicate numerical results. 

A = 0.00235, a = 12.1 nm, n = 1.457. The simplified theory shows good 
agreement. 

IV. CONCLUSION 

We have derived a very simple formula for the coupling between two 
multimode fibers and have discussed the coupling mechanism and the 
coupling efficiency. We emphasize that the formula obtained involves 
some rather drastic approximations. However, this coupling formula 
explains the coupling mechanism very clearly and agrees with experi- 
mental and numerical results. 
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APPENDIX A 
Coupling Coefficient 

Coupling between two modes of different fibers was discussed by 
Snyder. 1 - 2 We use the formula derived by Snyder. 
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We assume the two fibers to be nearly identical so that the pth mode 
of fiber A can only couple to the qth mode of fiber B, provided the two 
modes have the same z -direction propagation constant. 

In this case, the coupling equations have the following form 

dA(z) 



dz 



+ jf3 A A(z) = -jB(z)C AB 



^^ + j0 B B(z) = -jA(z)C B A, (15) 

dz 

where A (z) and B(z) are the amplitudes of the modes of fiber A and fiber 
B respectively. The coupling coefficient is defined by 

Cab -7; f <o(nl-n 2 c )E A E B dS, (16) 

Z %) core area 



where 



of fiber B 

no = core index of fiber B 

n c = cladding index of fibers A and B 

E A = the normalized electric field of fiber A 

E B = the normalized electric field of fiber B. 

Weakly guiding fibers have simple field expressions, which are called 
the linearly polarized modes (lp). Each LP[ m mode represents a set of 
four modes when / > 0. (When / = 0, LPo m represents two modes.) The 
four modes differ in polarization and azimuthal field distribution (the 
sin l(f> or cos l(f> term in the field-expansion equation). The field com- 
ponents can be described by 3 



E y - Hx 



Zq/uq 

zo/n c 



= A, 



Ji(ur/a)Ji(u) 
K t (wr/a)Ki(w) 



1 , . , ,0<rsa ,„„,. 

[cos 1(f) or sin t0j , (17) 

r > a 



where 



Al . /H) W (y/X) m^ I * 

W V v Ttno/ aivy/Ki-MKi+^w)] 

A„=(*) w (^)if^ ,-a 

MULTIMODE FIBER COUPLER 739 



These fields satisfy the following eigenvalue equation 
u[Ji-i(u)/Ji(u)] = -w[Ki- l (w)ki(w)) 



or 



u[Ji + i(u)/Ji(u)] = w[K l+l (W)/Ki(W)]. (18) 

We calculate the coupling coefficient between the pth mode of fiber A 
and the qth mode of fiber B, by using (16), (17) and (18). The pth mode 
field of fiber A is defined by h, w h a h r, 0, and the qth mode of fiber B 
is defined by l 2 , w 2l a 2 , r', 0'. The coupling coefficient is defined by (17). 

_ co € (n 2 -n 2 c ) 



J»2tt fCOS /i0 (COS/20')| 

o Isin h0 (sin l 2 <t>') \ 

x r w '*' i (<) j <°(" 2 3- <i9) 

We introduce the following theorem to change the coordinate (r 0) to 
the coordinate (r', 0') shown in Fig. 1: 

k„ U i) f cos ';* 1 = t <-i)*k» (». -) v* (». f) 

l \ ai/lsin/i0 a=— \ ai/ \ a x J 



The coupling coefficient is, 
co t (n|| - n c 2 ) 



icos[(/i + fe)0'-^a]. 
sin [(/1 + /z)0' - fea] 



p ' 2K|,(u)i)/(| 2 (w2) 

x J„ 02r ' dr ^(""3 j ' 2 ( u2 3 



2 X/ 1 (wi )K/ 2 (w 2 ) 

»?_„2 



(21) 



2K/,(«;i)K/ 2 (iy2) 



CO 6 (no - Mc)7T 

where 

Qtfi, h) = (-\) h ~ h K h -i 2 U>\—) cos (h - l 2 )a 

± (-l)ii+i2K ll+l2 ( Wl — ) cos (/1 + l 2 )a 
740 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1977 



/u 2 \ 2 /w 1 \2La2 \ o.\l 



©♦© 



+ —J[ 2 (u 2 )Il 2 + ! f iT, 



("'»} 



If we substitute the eigenvalue equation into R(li, l 2 ), we obtain 
a 2 Jl 2 {u 2 ) 



Rdul 2 ) = 



u 2\ 2 . /WA 2 K-hiwo) 



X 



©♦© 



— Ki 2+ i(w 2 )Ii 2 (wi—) 
a 2 \ a\l 

K)l- (22) 



+ — Ki 2 (w 2 )Ii 2+ i 



The final result for the coupling coefficient is 

kriQa 2 2A 2 u\U 2 



Ca p b, 



a\ V\v 2 



w 2 K[ 2+l (w 2 )Ii 2 + i (wi —J 



+ w l —Ki 2 {w 2 )Ii 2+ i (wi —) 
(-\)h-h Kh - h (wi -^j cos (/, - l 2 )a 

± (-l)'i+'2K/ 1+/2 ( Wl ^) cos (h + h)a 



(" 2 ) 2 +("i fj 



(K^-iiw^K^+iiwOK^+dw^K^-iiw^J 



1/2 



(23) 



APPENDIX B 

General Coupling Equation 

The coupling between two fibers, as shown in Fig. 1, is expressed by 
the following coupling equation with the coupling coefficient defined 
by (23). 



dA p (z) 
dz 



+ JP Ap A p (z) = -jj:B q (z)C ApBli 



^± + j0 Bq Bj(z) = -; Z A p (z)C BqAt 

CLZ n 
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( 7 i72) ^2AJ 

KW-— ^ ">" 2 



X 



Ki t -i 2 (wid/ai) cos (/i - l 2 )a ± gfr+jg (wid/ai) cos (ji + [g)a 



["■♦"*<Sfl 



(^2)^/2+ i(^2U/2(^i fl 2/qi) + wi(a2/ai)Ki 2 (wi)Ii 2+ i(wia 2 /ai) ( 2 4) 
[K Zl -i(«;i)K/ 1+ i(u;i)K, 2 _i(«; 2 )^ 2+ i(«;2)] 1/2 

where 

A p (z) = the amplitude of the pth mode of fiber A 

B q (z) = the amplitude of the qth mode of fiber B 

CApB q - coupling coefficient 

(3 A = the 2 -direction propagation constant of the pth mode in 

fiber A 
p B = the 2-direction propagation constant oftheqth mode in 

fiber B 
h = the azimuthal order number of the pth mode in fiber 

A 
l 2 = the azimuthal order number of the qth mode in fiber B 
Ki(z) = the modified Bessel function 
a 1 = the core radius of fiber A 
a 2 = the core radius of fiber B 
Ai = the normalized index difference between the core and the 

cladding of fiber A 
A 2 = the normalized index difference between the core and the 
cladding of fiber B 
71, 72 = when h = 0, l 2 = 0, 71 = 72 = 1; when h = 0,l 2 = 0, 71 = 
72 = 2. 
The power transfer from the pth mode of fiber A to qth mode of fiber 
B is given by the following equation where the pth mode of fiber A carries 
unit power and the qth mode of fiber B carries no power. 
\A p (z)\ 2 = 1 - k ab sin 2 PabZ 
\B q (z)\ 2 =K RA sin 2 Pabz, (25) 



where 



_r 1 , (fon-au 2 1- 1 

AB L 4|(C ApB J(C B , Ap )|J 



\ CbAp 

KRA ~ 



K-AB 



pAB _[ \(CA n B KC Bu .A D )\ V' 2 



I C Ap B q 

[ \(CA n B a )(CB u .A D )\ ~Y 
KAB J 

z = the coupling length of two fibers. 
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Strong coupling between two modes is expected only when the prop- 
agation constants are matched. Therefore the coupling between two 
multimode fibers, each of which carries hundreds of modes, can be an- 
alyzed by the coupling of modes whose propagation constants are 
matched. 



APPENDIX C 

Simplification of the Modified Bessel Function 

The asymptotic expansion for the modified Bessel function is given 
as follows 



tfifcj-Vf «-*i 



(l,n) 
2z~ „%(2z) n ' 



(26) 



where 



\z\ is large 
|argz | < 3tt/2 

(/ n) ~ (4/2 " 12)(4/2 " 32) ' * ' (4/2 - i2n " 1)2) _ V + " + 2/ 

»!r(l-„ + -) 

The parameters w and wd/a are large for most modes except those very 
close to the cutoff value. Therefore, most of the modified bessel functions 
can be expressed by the above asymptotic expansions. Especially when 
both order numbers of the mode are small (i.e., w is large), the modified 
Bessel function can be expressed by the first term of the asymptotic 
expansion. We obtain 



■ 2V 



K 2 i(wd/a) ±K (wd/a) 
_K^i(w)Ki- l (w) ' 

-— e-««/«-2) (« = 0) 

Trd/a 



s (a - 90°). (27) 

We show below that (27) can also be used for large /. The coupling 
coefficient has the following modified Bessel function term: 

K 21 (wd/a) ± K (wd/a) 

KmMKi-tiw) ' (28) 
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We combine (26) into (28); then the first term is, 

K 2l (wd/a) _ / 2w e _ w(d/a _ 2) 
Ki+i(w)K t -i(w) V nd/a 

- (2/, n) 



n % (2w) n (d/a) n 



- " (l+l,n-m)(l- 1, m) 

n =0m=0 (2W)» 



a \I ^L. e -w(d/a-2) X 



l„((i^)Q"[ 4(2 ' )2 ' n i 



ird/a °° 1 " 

£ tA: (4(/) 2 ) b E 



„to (2i») n m =o (n - m)\m\ 2 2 " 

f. J^(«\ n W) 2 \ n 22n 

+ I ^w ... nx n^o(2w)»\dJ 2 2 "n! 

= V e -w(d/a-2) m (29) 

V ird/a - 1 2"[4(/) 2 ]" 

n ?o(2^)" 2 2 "n! 

If d/a s 2, the above equation is 

K 2 i(wd/a) 
Ki + i(w)Ki-!(w) 

1 /a\" 2 2n (4l 2 ) n 



^ (2w) n W 



V 7rd/a 1 2"(4/ 2 ) n 

a */_?£' -«,<d/«-2) a (30) 
7rd/a 

Therefore, when / is large, the modified Bessel function term is expressed 
by the above equation. Then, 

K 2t (wd/a) ± Koiwd/a) _ -*/J^L e _ w(d / a - 2 ) (31) 

where 

£ = 1 when / large 

£ = 2 when / small, o> large (even) 

£ ^ when / small, u; large (odd). 
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Therefore, the coupling coefficient has the following limitation. 

2 5/4 A l/4 



0*\Ci\ iC B 



limit 



m-vr- - 



V irknoa 3/2 

However, we assume that all Bessel function terms are expressed by 
the case £ = 1 as the average coefficient. The result is 

.— , 2 1 / 4 A 1 ' 4 n 
where 



2 i / 4 a 1 / 4 rnr i v /4 
M-vj^-AnIV-n] • (33) 



C, = average coupling coefficient. 
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